Large Extra Dimensions, Sterile neutrinos and Solar Neutrino Data 
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Solar, atmospheric and LSND neutrino oscillation results require a light sterile neutrino, vb, which 
can exist in the bulk of extra dimensions. Solar v c , confined to the brane, can oscillate in the vacuum 
to the zero mode of vb and via successive MSW transitions to Kaluza-Klein states of vb ■ This new 
way to fit solar data is provided by both low and intermediate string scale models. From average 
rates seen in the three types of solar experiments, the Super-Kamiokande spectrum is predicted with 
73% probability, but dips characteristic of the 0.06 mm extra dimension should be seen in the SNO 
spectrum. 
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The four-neutrino scheme in which the solar v e deficit 
is explained by v e — > v s (where v s is a sterile neutrino), 
the atmospheric v^/v e anomaly is attributed to — ► v T , 
and the heavier and v T share the role of hot dark mat- 
ter was originally proposed H in order to explain those 
three phenomena. Later the LSND experiment , which 
observed — ► v e , provided a measure of the mass differ- 
ence between the nearly degenerate v e — v s and — v T 
pairs and required the three mass differences that were al- 
ready present in that neutrino scheme. Exactly this same 
pattern of neutrino masses and mixings appears neces- 
sary to allow production of heavy elements (A ~ 100) 
by type II supernovae While qualitatively this neu- 
trino scheme seems to explain all existing neutrino phe- 
nomena, solar neutrino observations are now sufficiently 
constraining that the small-angle MSW v e — > v s expla- 
nation appears to be in some difficulty Q. Although 
providing better fits to the solar data, even active-active 
transitions in a three-neutrino scheme do not give a quan- 
titatively good explanation of those data. In this Letter 
we point out that there is a way to achieve an excel- 
lent fit and rescue the apparently needed four-neutrino 
scheme using large extra dimensions. This is motivated 
by the latest developments in string theories, which have 
made plausible the interesting possibility that such large 
extra space dimensions (~mm) [^) can exist in conjunc- 
tion with low or intermediate fundamental string scales, 
M*. Using such models and the average rates from the 
three types of solar experiments, we predict the energy 
spectra from the Super-Kamiokande || and SNO experi- 
ments employing a combination of vacuum and MSW os- 
cillations. The neutrino-electron scattering results from 
Super-Kamiokande is predicted with a x 2 corresponding 
to 73% probability, while the better energy resolution of 
the charged-current interactions from SNO should show 
the predicted dips in the spectrum. 

We have found two classes of models giving the desired 
phenomenology (see jf)) for details). The first model has 
low M*, as has been used to give an alternative to the 



conventional high-scale SUSY GUT theories in solving 
the gauge hierarchy problem and providing low energy 
signals of string theories. Since there are no large scales 
to implement the conventional seesaw mechanism, a way 
to get small neutrino masses is to include singlet neutri- 
nos in the bulk, in combination with the assumption Rj 
of an effective global B — L symmetry in the theory below 
M* . The large size of the bulk suppresses neutrino masses 
to the desired level. The low Kaluza-Klein ("KK") exci- 
tations have a mass ~ ~ 10~ 3 eV, which not only 
provide a natural way to understand the lightness of the 
sterile neutrino ||, but also give extra sterile neutrino 
states in the mass range that can influence the shape of 
the solar neutrino spectrum ||. 

This model has one bulk neutrino, vs{x,y), y being 
the coordinate of the fifth dimension, and the standard 
model on the brand. The vb (x, y) is assumed to couple to 
the lepton doublet of the standard model L and of course 
have a five-dimensional kinetic energy term. After elec- 
troweak symmetry breaking at scale v w k, the v e — vb.r 



coupling leads to m v = h 



Mp t 



10~ 5 eV, where M Pe 



is the Planck mass and h is the Yukawa coupling. Note 
that this suppression is independent of the number and 
radius hierarchy of the extra dimensions, provided that 
ub propagates in the whole bulk. Even if the bulk is six 
or higher dimensional, there is only one mm-scale dimen- 
sion. The smaller dimensions will contribute only to the 
relationship between Mpi and M», but their KK excita- 
tions will be very heavy and decouple from the neutrino 
spectrum, making the analysis as in five dimensions. 

The direct Dirac or Majorana mass terms for the bulk 
neutrino can be forbidden by an appropriate choice of 
geometry and dimension of the bulk in which the vb 
resides, making an ultralight v s natural. For instance, 
in 5 dimensions the orbifold symmetry under which 
y — > — y combined with B — L symmetry guarantee this 
for the conventional definition of charge conjugation. 

In order to fit neutrino data, one needs to include new 
physics in the brane that will generate a Majorana mass 
matrix for the three standard model neutrinos of the form 
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5 a b (where a, b = e, p, r). For 5 MT much bigger than the 
other elements, the ^ ljT in effect decouple from the v etS 
and do not affect the mixing between the bulk neutrino 
modes and the v e . Further, this leads to maximal mixing 
in the fi — r sector, as is needed to understand the at- 
mospheric neutrino data. If 8^ T ^eV, then this provides 
an explanation of the LSND observations. One way to 
generate this pattern is to consider an L e +L^ — L T sym- 
metric extension of the standard model with additional 
doubly-charged singlet (Y = 4) and Y — 2 triplet scalar 
fields. In the rest of the Letter, we focus only on the 
v e — v s sector and how we fit the solar neutrino data. 
The mass matrix for the v e , v s sector can be written 
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where v' B represents the KK excitations, and the off- 
diagonal term v2to is actually an infinite row vector of 
the form ^2^(1, 1, ■ • ■). The operator d$ stands for the 
diagonal and infinite KK mass matrix whose n-th entry 
is given by n/R. Using this short-hand notation makes 
it easier to calculate the exact eigenvalues m n and the 
eigenstates of this mass matrix, 



m n = S ee + nm 2 Rcot(Trm n R). 
The equation for eigenstates is 

a/2to {m n v' B _ + d 5 is' B + ) 
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where the sum over the KK modes in the last term is 
implicit. N n is the normalization factor given by 
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Note that in the limit of S ee = 0, the v e and pb,o are two, 
two-component spinors that form a Dirac fermion with 
mass m. The KK modes come in pairs of mass m n = 
±fc„/i?, with k n a positive integer, and they couple to 
the v e approximately as mR. Once we include the effect 
of See 7^ 0, they become Majorana fermions with masses 
given by to-i ~ +<5 ee /2 + to and to 2 ~ +S ee /2 — to, and 
they are maximally mixed; i.e., the two mass eigenstates 
are v\.2 — V '^j^' a ■ Thus as the v e produced in a weak 
interaction process evolves, it oscillates to the vb,o state 
with an oscillation length D ~ E /(2mS e e), which for 
natural values of m, S ee gives D of order of the Sun-Earth 
distance so that our model leads to vacuum oscillation 
("VO") of the solar neutrinos. Furthermore, since the v e 
also mixes with the KK modes of the bulk neutrinos with 
a 8m 2 ~ 10 -5 eV 2 , this brings in the MSW resonance 
transition of v e to vb^kk modes at higher energies. 

The second way to achieve the same phenomenology 
is to use a much higher string scale associated with local 



SU(2) L x SU(2) R x U{\) B -l symmetry of the left-right 
symmetric model in the brane, coupled to a single vb in 
the large fifth dimensional bulk. Using the usual fermion 
content of the left-right models and breaking the right- 
handed symmetry by Higgs doublets, one obtains the see- 
saw type matrix in the presence of the infinite tower of 
KK modes. The heavy right-handed neutrino now de- 
couples, leading to the following neutrino mass matrix in 
the one generation case (the basis is {v e VQBv' B -v' B +)) 
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(1) where a 



hi_M*VR 



= hpv wk , M = v R /M*, and v R is 
the scale of SU{2) R breaking. A denotes the radiative 
correction induced due to the extrapolation from the v R 
scale to the weak scale. The lowest eigenvalue of this 

2 A 

mass matrix is ~ ../"f and the next lowest one 

is m j +" . For to, a ~ 1 — 5 MeV (similar to the first 
generation fermion mass) and M ~ 10 9 GeV, we get this 
eigenvalue to be of order 1 0~ 6 — 2.5 x 10 -5 eV. Its square 
is therefore in the range where the VO solution to the 
solar neutrino puzzle can be applied. Also for to ~ a, 
the mixing angle between the zero eigenvalue mode and 
this mode is maximal. Thus this model has properties 
similar to the first model for neutrinos, and below we 
carry out our fit to solar neutrino data using the latter. 

For propagation in solar matter, the eigenvectors and 
eigenvalues can be found by replacing the squared mass 
matrix, M 2 , for the neutrinos, with M 2 + H, where H = 
2Ep e when acting on v & (where p e = \f2Gpine — n n /2)), 
and zero on sterile neutrinos. Defining 
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(w = 1 in vacuum) Eq. || becomes 

m n — wS ee + irm 2 R cot(nm n R) . 



(7) 



The eigenvectors are as in Eq. |[ except the third term 
acquires an additional factor l/w. 

In fitting the solar data using VO, the strategy gen- 
erally employed is to suppress Be neutrinos while also 
reducing the 8 B neutrinos by more than 50%. The wa- 
ter data then requires an additional contribution, which, 
in the case of active VO, is provided by the neutral cur- 
rent cross section amounting to about 16% of the charged 
current one. Thus in a pure two-neutrino oscillation pic- 
ture, VO does not work for active to sterile oscillation. In 
our model, however, both vacuum oscillations and MSW 
oscillations are important, since the lowest mass pair of 
neutrinos is split by a very small mass difference, whereas 
the KK states have to be separated by > ICP 3 eV because 
of the limits from gravity experiments. We can use the 



first node of the survival probability, P ee , to suppress 
the 7 Be. Going up in energy toward 8 B neutrinos, P ee , 
which in the VO case would have risen to very near one, 
is reduced by the small-angle MSW transitions to the 
different KK excitations of the bulk sterile neutrinos, as 
is clear from Fig. [I]. This is a new way to fit the solar 
neutrino data in models with large extra dimensions. 

To do the fit, we studied the time evolution of the v e 
state with a program supplied by W. Haxton G3], but up- 
dated to use a recent solar model Jll|] and modified to do 
all neutrino transport within the Sun numerically, using 
no adiabatic approximation. Changes were also neces- 
sary for oscillations into sterile neutrinos and to general- 
ize beyond the two-neutrino model, for up to 14 neutrinos 
contribute for the solutions we considered. 
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is 0.451 ± 0.016. The best fits were with R « 58 fim, mR 
around 0.0094, and S ee ~ 0.84 x 10~ 7 eV, correspond- 
ing to Sm 2 ~ 0.53 x 10~ n eV 2 , giving average P ee for 
Chlorine, Gallium, and water of 0.383, 0.533, and 0.450, 
respectively, and the P ee energy dependence shown in 
Fig. fy. For two-neutrino oscillations, the mixing angle 
is sin 26, whereas here the coupling between v e and the 
first KK excitation replaces sin 2 28 by Am 2 R 2 — 0.00035. 

Vacuum oscillations between the lowest two mass 
eigenstates nearly eliminate electron neutrinos with en- 
ergies of 0.63 MeV/(2n+l) for n = 0, 1, 2, ... . Thus 
Fig. [l] shows nearly zero P ee near 0.63 MeV, partly elim- 
inating the 7 Be contribution at 0.862 MeV, and giving a 
dip at the lowest neutrino energy. MSW resonances with 
mass pairs of higher KK states start causing the third 
and fourth eigenstates to be significantly occupied above 
~ 0.8 MeV, the fifth and sixth eigenstates above ~ 3.7 
MeV, the 7'th and 8'th above ~ 8.6 MeV, and the 9'th 
and 10'th above ~ 15.2 MeV. Fig. |l] shows dips in P ee 
just above these energy thresholds. 

The expected energy dependence of P ee is compared 
with Super-K data [Q in Fig. ^. The uncertainties are 
statistical only. The parameters used in making Fig. |^ 
were chosen to provide a good fit to the total rates only; 
they were not adjusted to fit this spectrum. Combining 
spectrum data with rates using the method described 
in Ref. jl6) gives \ 2 — 14.0 for the spectrum predicted 
from the fit to total rates. With 18 degrees of free- 
dom, the probability of x 2 > 14.0 is 73%. A fit with 
5 ee constrained to be very small to eliminate vacuum os- 
cillations increased the best fit \ 2 from 3.4 to 4.4. The 
same parameters then used with the Super-K spectrum 
gave x 2 = 18-7 (probability 41%). This is comparable to 
\ 2 = 19.0 for an energy independent spectrum. 



FIG. 1. Energy dependence of the v e survival probability 
when R = 58 fj,m, mR = 0.0094, and 5 ee = 0.84 x 10" 7 eV. 
The dot-dashed part of the curve assumes the radial depen- 
dence in the Sun for neutrinos from the pp reaction, the solid 
part assumes ls O radial dependence, and the dashed part as- 
sumes S B radial dependence. 

For comparison with experimental results, tables of de- 
tector sensitivity for the Chlorine and Gallium exper- 
iments were taken from Bahcall's web site |TlJ. The 
Super-Kamiokande detector sensitivity was modeled us- 
ing for the resolution in the signal from Cherenkov 
light. More details are given elsewhere ||. Calculations 
of P ee , averaged over the response of detectors, were com- 
pared with measurements. While theoretical uncertain- 
ties in the solar model and detector response were in- 
cluded in the computation of \ 2 as described in Ref. 
fL3|| , the measurement results given here include only 
experimental statistical and systematic errors added in 
quadrature. The Chlorine P ee from Homestake is 
0.332 ± 0.030. Gallium results fl§ for SAGE, GALLEX 
and GNO were combined to give a P ee of 0.579 ± 0.039. 
The 5.0-20 MeV, 1258-day Super-K experimental P ee |] 
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FIG. 2. Super-Kamiokande 1258-day measured j|] energy 
spectrum (error bars) and predicted (curve) using the param- 
eters of Fig. |l| without fitting these data. 
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The seasonal effect was computed for a few points on 
the Earth's orbit. If r is the distance between the Earth 
and the Sun, ^ = 1 + ecos(9 — 8 ), where r is one 
astronomical unit, e = 0.0167 is the orbital eccentricity, 
and 9 — 0Q ~ 2ir(t — to) , with t in years and to = January 
2, 4h 52m. Table I shows very small seasonal variation. 

Not only is the seasonal effect very difficult to observe, 
but also the smallness of the mixing angle makes day- 
night effects hard to measure. In addition, the mass 
of the electron neutrino, which consists mainly of eigen- 
states of mass 3 x 10~ 5 eV, is undetectable directly or by 
neutrinoless double beta decay. The latter process mea- 
sures an effective neutrino mass, but even contributions 
to that from the and v r must be so small as to make 
detection very unlikely, although other conjectured pro- 
cesses unrelated to neutrino mass could cause this decay. 

On the other hand, the dimension size of 0.06 mm, 
suggested by the average rates of the three types of solar 
experiments should be detectable by gravity experiments. 
The present best limit |l7j on such effects is less than a 
factor of four from that value. 

The large size of the extra dimension raises issues 
about cosmological and supernova limits from the effects 
of high KK states of both sterile neutrinos ]l^-[20|] and 
gravitons pjj , despite uncertainties in the understand- 
ing of the complex regimes of the early universe and 
the supernova core. For sterile neutrino limits, this phe- 
nomenology is aided because there is a single KK tower 
based on an exceedingly small mass, the VO Am 2 is an 
order of magnitude smaller than usual, and for MSW 
the equivalent sin 2 29 is more than an order of magnitude 
smaller than for standard fits. For the global B-L model, 
the universe re-heat temperature could be very low (> 0.7 
MeV works cosmologically), reducing production of high 
KK states. The high string scale of the local B-L model 
would appear to avoid all these constraints how- 
ever. More complete investigation of these constraints 
may enable choosing between these quite different mod- 
els, both of which provide this new way to explain solar, 
atmospheric, and LSND oscillation data and may give 
the first evidence for an extra large dimension. 
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